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ABSTRACT 

We analyse new integral field spectroscopy of the inner region (central 2.5 kpc) of the 
spiral galaxy NGC 4321 to study the peculiar kinematics of this region. Fourier anal- 
ysis of the velocity residuals obtained by subtracting an axisymmetric rotation model 
from the Ha velocity field, indicates that the distortions are global features generated 
by an m = 2 perturbation of the gravitational potential which can be explained by the 
nuclear bar. This bar has been previously observed in the near-infrared but not in the 
optical continuum dominated by star formation. We detect the optical counterpart of 
this bar in the 2D distribution of the old stellar population (inferred from the equiv- 
alent width map of the stellar absorption lines). We apply the Tremaine- Weinberg 
method to the stellar velocity field to calculate the pattern speed of the inner bar, 
obtaining a value of r2i,=160±70 km s _1 kpc . This value is considerably larger than 
the one obtained when a simple bar model is considered. However the uncertainties in 
the pattern speed determination prevent us to give support to alternative scenarios. 

Key words: galaxies: individual (NGC 4321) — galaxies: kinematics and dynamics. 



1 INTRODUCTION 

Departures from regular rotation, such as winds, elongated 
orbits or streaming motions, provide very relevant informa- 
tion about the dynamics and history of galaxies. To identify 
kinematical distortions and to infer how they relate to the 
gravitational potential, different techniques such as velocity 
field fitting or Fourier analysis can be used. However, to ob- 
tain reliable results these tools must be used in combination 
with 2D kinematical data of good spatial resolution. Inte- 
gral field spectroscopy (IFS) is an observational technique 
especially suited for kinematical studies that can be applied 
to nearby galaxies, such as M100 (NGC 4321), to obtain 
velocity maps with detailed spatial information. 

The inner region of NGC 4321 has been exten sively 
studied at di fferent wavelengths. O ptical (|Pierceil 19861). Ha. 
near- in frared ([Knapen et al.lll 995a b) a nd CO studies |Rand| 

irilld 



near-infrared shows an inner bar aligned with the large 
scale ste llar bar and a pair o f small arms emerging from 
its ends (|Knapen et al.|[l995al fbh. Whether there is a single 
bar or two nested bars that just happen to be aligned 
has b een a great l y debated issue |Knapen et alj Il995bl . 
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1995 : ISakamoto et al.l 1 19951 : Isempere fc Garcfa-Burillc 
19971 ). have revealed some interesting irregular features. 
On the one hand, the inner region of this galaxy shows a 
peculiar morphology that is quite different in the optical 
and the NIR. In the optical, the structure is dominated 
by two spiral arms (broad band) and an ovally shaped 
region of enhanced star formation (Ha, H/3), while the 
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data, and 

kinematical signatures of gas streaming along the inner 
part of the bar, and of density-wave streaming motions 
across a two-armed minispiral. Relevant pattern speeds of 
this gala xy have been measured at larger scale by different 
meth o ds (ICanzian fc Allenlll997l ; ISempere fc Garcfa-Burillol 
1 19951 : iHernandez et al.| [20051. _ and references therein). 
Recently, lAllard et al.1 l|2005al ) have analysed the inner 
region of NGC 4321 from IFS observations obtained with 
SAURON. They focus on the star formation ring, which 
they confirm to be formed by young stars, but find a low 
gas velocity dispersion compared to its surroundings. 

In this paper we present new IFS (covering a large wave- 
length range in the optical) obtained with INTEGRAL, a 
system of exchangin g fibre bundles at the WHT suited for 
this kind of study dMediavilla et al.1 [2005; Bat taner et al.1 
120031 : iGartia-Lorenzo et al.1 1200 li b We use the stellar and 
ionized gas velocity fields to analyse the symmetries of the 
departures from pure rotation and to estimate the nuclear 
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Name 


M100 


Type 


SAB(s)bc LINER-HII 


R.A. (J2000) 


12 fe 22 m 54.9 s 


Dec. (J2000) 


15°49'20" 


Incl. 


30° 


P.A. 


151° ±3° 


Total Mag. (B) 


10.05 


Sys. Velocity (opt) 


1567 ± 7 km/s 


Distance 


16.1 Mpc 


Linear scale 


78 pc/arcsec 



Table 1. M100 main properties 



bar pattern speed by means of the Tremaine & Weinberg 
(1984) method. These studies will be used to relate the mor- 
phological and kinematical distortions observed in the inner 
2.5 kpc with the gravitational potential of the nuclear bar. 

NGC 4321 has a very accurately (Ceph eid) measured 
distance of 16.1 Mpc l|Ferrarese et alj Il996l ) which corre- 
sponds to 78 pc arcsec -1 on the plane on the sky. We use 
this value throughout this work. 

The main properties of M100 are summarized in Ta- 
bled] 



2 OBSERVATIONS AND DATA REDUCTION 

The data analysed in this article were obtained on 2002 
March 16 at the Observatorio del Roque de los Mucha- 
chos on th e island of La Palm a with the fibre system IN- 
TEGRAL (|Arribas et al.lll998l ) in combination with the fi- 
bre sp ectrograph WYFFOS at the William Herschel Tele- 
scope (iBingham et al.lll994l ). The weather conditions dur- 
ing this night were fairly good, with a seeing of about 1.3". 
The data discussed in this paper were obtained with INTE- 
GRAL standard bundles #3 and #2. The WYFFOS spec- 
trograph was equipped with a 1200 groove mm -1 grating 
centred on 6247 A. The spectral resolution was 4.8A (R w 
1300) for SB3 and 2.8A (R w 2200) for SB2 fibre bundles. 
INTEGRAL-I- WYFFOS allow us to observe a large (com- 
pared with other IFS systems) spectral range (5600-6850A) 
that, in addition to Ha, also contains the high excitation [N 
IIJAA6548, 6584 and [S IIJAA6716, 6731, emission lines and 
the Na D AA5890, 5896 absorption features (see Fig. ED (a)) 
not included in previous IFS studies. With this configura- 
tion, and pointing to the centre of NGC 4321, we took three 
exposures of 1200 s each during the night. 

The data reduction of fibre-based 2D spe ctroscopy ( 
e.g. lArribas et ail Il99ll ; iMediavilla et all 1 19921 ) consists of 
two main steps: i) basic reduction of the spectra (i.e. bias, 
flat-fielding, extraction, wavelength calibration, etc.) and ii) 
generation of maps of spectral features (e.g. line intensity, 
velocity fields, etc.) from the spectra. Step (i) was performed 
in the IRAF environment following standard procedures. 
We obtained typical wavelength calibration errors of 0.15 
A, which give velocity uncertainties of ± 7 km s _1 for Ha. 
For step (ii) we have developed our own software packages. 
In particular we transform an ASCII file with the actual po- 
sition of the fibres and the spectral feature corresponding to 
each fibre into a regularly spaced rectangular grid. In this 
way we build up images of 33 x 29 pixels for SB3 bundle, 
and 45 x 34 pixels for SB2, with a scale of ~ 0.95" pixel -1 , 



and 0.35"pixel for SB3 and SB2 fibre bundles respectively. 
The resulting images are therefore of sizes 31.35" x 27.55" 
for SB3 and 15.75" x 11.9" for SB2. These images can be 
treated with standard astronomical software. 



3 RESULTS 

3.1 Nuclear Spectrum 

The stellar population in t he central region o f NGC 4321 has 
been carefully analysed bv lSarzi et alj ([2005). They estimate 
an average age of 1 Gyr at solar or super-solar metallici- 
ties, with only a marginal improvement in the results for 
super-solar metallicities. We therefore use a synthetic stel- 
lar population of 1 Gyr with solar metallicity as a template 
for the stellar absorptio n feature analys i s. We have used the 
synthetic population bv;Del gado et aL~l |2005l) . which covers 
the observed wavelength range at a very high spectral reso- 
lution (0.3 A). This range includes the strong NaD line that 
might be affected by interstellar absorption. We thus exclude 
the wavelength range of NaD from the fitting, obtaining a 
template that reproduces very well the observed absorption 
features, except the NaD line, which exhibits strong inter- 
stellar contamination. The difference between the nuclear 
spectrum and the template shows that the interstellar con- 
tribution is blue-shifted (see Fig. [Tj (b)). This is a notable 
result, since blue-shifted velocity components in this fea- 
ture unambiguously indicate the presenc e of outflowing gas , 
as demonstrated in previous studies (see iRupke et al.ll2005l . 
and references therein). 



3.2 Ha Maps 

To obtain the intensity and velocity maps for the ionized gas 
we fit a single gaussian to the Ha emission line after sub- 
tracting the stellar contribution from each spectrum. Ab- 
sorption correction is in general not very important except 
at the nucleus. According to the Ha intensity map (see Fig. [2] 
(a)), the em ission is distributed in a nucleus and a n ovally 
shaped ring l|Knapen et al ] |2000l ; lAllard et al . 2005a), where 
several HII regions can be distinguished. 

The Ha velocity field map is computed with the com- 
bined information from the SB2 and SB3 fibre bundles (see 
Fig[3](a)). O ur velocity map looks quite similar to the one 
published by iKnapen et al.l (|2000l ). It essentially shows a 
basic rotational pattern significantly distorted by the effect 
of the inner bar and spiral arm structures. In FigO (b) we 
overplot the ionized gas velocity contours on the HST/ACS 
(F555W) image. In this figure, the two regions with clear 
departures from axisymmetric circular rotation stand out 
clearly. The innermost region shows a clear twist in the kine- 
matic axes due to the influence of the inner bar. Further out, 
the deviations follow the spiral arms. 



3.3 Stellar Maps 

Cross-correlation of the spectra is performed in the 5600- 
6560A wavelength range with the XCSAO task from the 
IRAF package. We have excluded in the cross-correlation 
wavelength ranges around the [O I]AA6300, 6363 lines to 
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Figure 1. (a) Nuclear spectrum (fibre #68 in bundle SB3) of NGC 4321 in the range 5600—6850 A showing the most important 
absorption and emission lines, (b) Nuclear spectrum for the SB3 fibre bundle with the stellar template plotted on top in red line. This 
template is used with the XCSAO task to derive the stellar velocity field using the cross-correlation technique in the wavelength range 
5600-6560A. The inner panel shows in detail the fit of the Na D absorption doublet (red line) and the residuals (blue line). 
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Figure 2. (a) Ha intensity map computed for the SB3 fibre bundle. The Ha emission is distributed in a nucleus and an ovally shaped 
ring, (b) HST/ACS (F555W) image of NGC 4321 with white contours overlaid from the continuum map computed for the SB2 fibre 
bundle. 



avoid residual sky emission contamination and a wave- 
length range around the NaD feature (for the reasons men- 
tioned above) . The stellar velocity uncertainties ranges from 
15 km s _1 (inner regions) to 35 km s _1 (outer regions). The 
resulting stellar velocity field is shown in Fig. [4] (a) . It look s 
very similar to the results published bv lAllard et alJ (|2005bl ) . 
Although noisier than the gas velocity field, it shows a stel- 
lar rotation pattern that is more regular than the velocity 
pattern of the ionized gas. It shows some distortions in the 
bar region although not as strong as in the case of the gas. 
The strongest deviations from circular rotation in the map 
take place at two spots in the spiral arms. They are located 
roughly symmetrically with respect to the galactic centre 
with a P.A. of about 45° from the major axis (nearly perpen- 
dicular to the bar P.A.), c oincident with the high est density 
regions of the CO map bv lSakamoto et all (|l995t) . This sup- 



ports the role of the nuclear bar in explaining the complex 
kinematics and morphology in the innermost region of NGC 
4321. This issue, and the question of whether the nuclear 
bar is part of the large scale one present in th is galaxy or a 
separated structure (see lHernandez et al.ll2005l . for a review) 
will be addressed in Section [ 



The stellar continuum map shows very similar features 
to the broad band HST/ACS (F555W) image although at a 
lower spatial resolution (see Fig. H](b)). In this map we can 
distinguish the nucleus emission from which a spiral struc- 
ture is coming out in the SW direction and a fainter one at 
the NE direction. Dust structures at the inner rims of the 
arms are als o clearly seen in the INTEGRAL continuum 
map (see also lKnapen et al1ll995al lbl). The brightest regions 
in this map trace stellar formation regions, of which is par- 
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Figure 3. (a) Ha velocity map computed with the information from the SB2 and SB3 fibre bundles. Contour levels range from 1470 
km s — 1 to 1650 km s — 1 in steps of 10 km s — 1 . The thick white contour indicates the systemic velocity at 1567 km s . Figures (a) and 
(c) share the same scale colour bar. (b) HST/ACS (F555W) image of NGC 4321 with white contours of Ha velocity field overlaid, (c) 
Model velocity map with fixed PA (see text), (d) Map of residual velocities derived by subtracting the model velocity field (c) from the 
observed velocity map (a). Contour levels: from —30 km s — 1 to 50 km s — 1 in steps of 8 km s — 1 . Yellow and red lines show the position 
angle of the galaxy (P.A = 151°) and the inner stellar bar (P.A = 111°) respectively. 



ticularly clear the bright region surrounded by dust in the 
NW of the nucleus. 

We have also calculated maps for the strength and 
equivalent width of the stellar absorption lines from the fit- 
tings of the template to the spectra obtained during the 
cross-correlation procedure (see previous section). It should 
be kept in mind that we have used the same single age, sin- 
gle metallicity population template for the whole map, and 
that therefore the maps for the intensities and equivalent 
widths of all the stellar absorption lines are proportional. 
This is obviously an oversimplification, as we know there 
will be variations in stellar populations throughout the field 
of view, but a thorough analysis of the stellar populations is 
beyond the scope of this paper, and these maps are intended 
only to be a tracer of the global stellar mass content of the 
galaxy. Notice that we have obtained a map normalized (ac- 



cording to the template) to the NaD line absorptions, which 
are the strongest ones in the wavelength range studied. This 
NaD map is therefore free from interstellar contamination. 

Figure [5] shows the contours of the equivalent width 
of t he stellar NaD line on top of the WHT Ks image 
from iKnapen et all (|2003l ). The iso-levels of the EW in the 
innermost region a re elongated with a P.A. very close to 
the 111° quoted bv lKnapen et all ljl995at ) for the inner bar. 
The extent of the bar does also r oughly agree with the 9 
arcsec proposed bv lKnapen et all (|l995al ). Although this is 
not completely surprising (as the equivalent width is ex- 
pected to be somewhat less affected by dust than contin- 
uum or broadband maps), it clearly shows that the EW is 
a good tracer of the old population responsible for the ab- 
sorp tion features. The two well known infrared knots Kl and 
K2 (jKnapen et al.|[l995al ) do not show up in the EW map 
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Figure 4. (a) Stellar velocity field for the SB3 fiber bundle. Velocities range from 1470 km s —1 to 1650 km s — 1 . The thick black contour 
indicates the systemic velocity at 1567 km s _1 . Yellow and red lines show the position angle of the galaxy (P. A. = 151°) and the inner 
stellar bar (P.A. = 111°) respectively. Contours of the CO emission are overplotted in white, (b) Stellar (star symbols) and Ha (diamond 
symbols) rotation curves of the inner disk of NGC 4321 derived up to 12" (fa 1 kpc) from the galactic centre. 




12:22:55.0 
R. A. (J2000} 



Figure 5. Infrared (K s ) image of NGC 4321 from jKnapen et al.l 
(2003J) with contours of the NaD equivalent width for the SB3 
fibre bundle. Contour levels range from 0.75 to 1.3A in steps o f 
0.05A. Infrared knots Kl and K2 dKnapen et aT I ll995ah 
are marked. 



and in fact correspond to minima on it (see Figure [5|). 
These knots a re known to be regions of powerful recent 
star formation IjRvder fc Knapenll999l ; IWozniak et al.lll998T ) 
and are therefore very bright because of the contribution of 
young stars, but not necessarily rich in old stars respon- 
sible for the absorption features contributing to the EW. 
Although the spectra are smeared by poor angular resolu- 
tion, we have been able to detect, for the first time, the 
(otherwise hidden) inner bar through the study of absorp- 
tion spectral features characteristic of old stars. We would 
like to point out that, t he K s band image sho ws very little 
hint of the spiral arms (|Knapen et al.lll995al lbl). This fact 



indicates that the mass distribution is quite axisymmetric 
at radii larger than « 6". 



4 DISCUSSION 

4.1 Residual Velocity Field of the Ionized Gas 

Departures from circular motion are better addressed 
through residual velocity analysis. We therefore start by de- 
riving a pure rotation model by fittin g tilted rings to the 
observed velocity field (|Begemanl ll989h. This is a standard 
procedure for which the galaxy is divided into concentric 
elliptical annuli where each annulus is characterized by the 
inclination angle (i), the position angle of the major axis 
(P.A.), the rotational velocity (v c ), the systemic velocity 
of the galaxy (v sya ) and the coordinates of the centre of 
each annulus (x c ,y c ). These parameters are fitted following 
a least squares algorithm. We have used the task ROTCUR 
from the data reduction software package GIPSY (Gronin- 
gen Image Processing SYstem) for this purpose. 

We take the kinematical centre at the position of the nu- 
clear source at our highest available spatial resolution. This 
position matches, within the resolution limits, the position 
of the brightest spot in the HST image. We fix the inclina- 
tion to a value of 30° dKnapen et al.l [20021. The systemic 
velocity is calculated as the average value of the outermost 
rings. We obtain a value of 1567 ± 7 km s _1 . With these 
values fixed, we perform a fit leaving the P.A. and the rota- 
tional velocity of each ring as free parameters. The results 
are shown in Figure [U The error bars in the rotation curves 
of Figures [U] and [3] b) are calculated from the difference 
in the rotation curves calculated for the approaching and 
receding sides of the galaxy and the error in the velocity 
determination, which have been added quadratically. 

The rotation curve obtained s hows a very similar 
behaviour to previous observations (|Knapen et al.l |2000| ; 
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Figure 6. Ha rotation curve of the inner disc of NGC 4321 de- 
rived up to 12" (f» 1 kpc) from the galactic centre. The lower 
panel shows variations of P.A. with radius. 



ISakamoto et al. 1 1 19951 ). Our values are in fact somewhat 
lower than those reported by iKnapen et al. | l|2000t). but in 
pod agreement with CO observations by ISakamoto et alj 
19951 ). Although, as usual, the inner region provides the 
most uncertain determination of the rotation curve, we have 
checked that rotations curves calculated for the maps at the 
two available resolutions (SB2 and SB3 fibre bundles) give 
an excellent agreement for the values of the rotational veloc- 
ity in the overlapping regions. Resolution effects will there- 
fore not be too relevant in the shape of the inner rotation 
curve. 

The rotation curve shows a very steep rise in the inner 
3", reaching values of 135 km s _1 . The velocity reaches a 
maximum of about 160 km s _1 at 6" (470 pc) and then 
declines smoothly to a value of 130 km s _1 in the outermost 
ring at 12" (« 1 kpc). 

To investigate the non-circular motions in the galaxy 
disc we compute the residual velocity map. We would like to 
minimize the effect of non-circular motions being absorbed 
in the rotational model. We use for this purpose a rotational 
model with a fixed P.A. (= 151° ± 3°, which is the average 
of the P.A. values corresponding to the outermost rings). 
In this way, most of the information on non-axisymmetric 
components is kept in the residual velocity. Of course, devi- 
ations from circular rotation that take place along or near 
the major axis in a symmetric way (on both the approach- 
ing and the receding sides) will most probably be absorbed 
in the rotation curve shape. There is little hope of disen- 
tangling the non-circular and circular components along the 
major axis, and therefore we adopt this as the best possible 
solution. The calculated axisymmetric rotational model is 
displayed in Figure [3] (c). Subtracting this model from the 
observed velocity field, we obtain the residual velocity map 
for NGC 4321 which is shown in Figure |3] (d) . 

Three regions can be clearly identified in the residual 
velocity map: (i) A blue-shifted nuclear region that shows 
blue residuals of about —25 km s _1 . We would like to point 
out that although the rotation curve is quite uncertain in 
the innermost region, the existence of these blue residuals 



in the nuclear region does mostly depend on the position 
of the kinematic centre and on the systemic velocity, which 
are both pretty well constrained by the observations. This 
fact suggests the existence of an outflow in this region, (ii) 
A second region, further out up to radii of about 6", that 
shows a clear signature of influence of the nuclear bar (whose 
position angle and size are shown in Fig. [3] (d) by the red 
line). The map shows residuals of nearly 50 km s _1 on the 
west side and of —25 km s _1 on the east side. These are 
the zones where orbits tend to align with the bar potential. 
The symmetric shape of this region above and below the 
bar is quite remarkable, (iii) At larger radii (r > 6") the 
residual map shows clear signatures of the presence of the 
spiral arms with values of 30 km s _1 for the northern arm 
and —25 km s" 1 for the southern arm. The northern arm 
is more clearly depicted than the southern one. In fact, as 
can be seen in Fig. [7] (a), the residuals for the northern 
arm nicely mat ch the position o f the molecular gas from 
the CO map by ISakamoto et all |l995). The southern arm 
is somewhat displaced from the position where the highest 
residuals are found. The residual velocities along the arms 
seem to vanish as they approach the major axis, although, 
as we have mentioned before, it is very likely that a large 
fraction of the non-circular motions has been absorbed in 
the shape of the rotation curve. The l atter two of the thre e 
regions have already been identified bv lKnapen et al.l {2000). 
In Section [4.21 we perform a Fourier analysis of the ionized 
gas kinematics to study the relationship of these regions with 
the gravitational potential. 

It is worth noting that the highest Ha intensity spots 
take place in reg ions of low residual velocities (see also 
IZurita et ai1l2004r i. In fact, in the case of NGC 4321, the Ha 
ring of star formation actually lies in the transition region 
between the bar and arm regions (with a radius of about 
6") with very low residual velocities. It is indeed very likely 
that both facts (high Ha intensity and low residual veloci- 
ties) are closely related and star formation takes place in a 
kinematically quiet region. 

A careful look at these residuals also reveals that the 
highest velocity gradients (which trace shocked gas regions) 
take place at the inner rim of the spiral arms and around 
the nucleus of the galaxy (at its far side). This is better 
illustrated if we calculate the gradie nt of the resid ual veloc- 
ity field. We have used the Roberts (|Robertsll 19651 ) operator 
to compute the 2-D spatial gradient measurement on the 
residual velocity field. The residual velocity gradient con- 
tours are shown on top of an unsharp-masked HST image 
of NGC 4321 in Figure [7] (b). It is quite impressive how 
the shocked gas traced by the high gradient matches very 
well the location of the dust lanes along the inner rim of the 
arms (particularly for the northern arm). This indicates that 
the spiral arms are very probably formed by hydrodynamic 
shocks originated by gas orbit crowding at a resonance of 
the bar potential (as the mass distribution as traced by the 
NIR light shows no relevant overdensities at these locations) . 
This coincidence of the dust location with high velocity gra- 
di ents was also found and studied in detail for NGC 1530 
bv lZurita etafl (|2004l ). 

The nuclear velocity gradient, although mixed with 
some bar-generated shocks, may also be understood in the 
context of the presence of an outflow in the nuclear region. 
If the approaching side of an outflow (more likely to be de- 
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Figure 7. (a) The residual velocities (colour image) nicely match the position of the molecular gas from the CO emission shown in 
white contours, (b) HST/ACS (F555W) unsharp masked image of NGC 4321 with residual velocity gradient contours overlaid. Yellow 
and red lines show the position angle of the galaxy (P.A. = 151°) and the inner stellar bar (P.A. = 111°) respectively. 
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Figure 8. (a) Residual velocity field (colour image) with the symmetric and antisymmetric dominance regions overplotted on white and 
black contours respectively, (b) Radial variation in amplitude of the Fourier modes of the residual velocity field. We show the results for 
the first four modes. 



tected than the receding one, located behind the galactic 
disc) produces a region of shocked gas over the nuclear re- 
gion, this would be seen in projection somewhat displaced 
along the minor axis on the far side of the disc. Additional 
indications of an outflow can be found in the morphology 
of the emission line profiles of the innermost region spec- 
tra that show blue shoulders and blue wings, especially 
in the h igh excitation lines ([Nil], S HI). In other galaxies 
(see e .g. lArribas fc Mediavillal (|l99a ), lAr"ribas fc Mediavillal 
(1994)). these features are an evidence of streaming gas. Fi- 
nally, the existence of a nuclear outflow is also supported 
by the blue-shifted interstellar contamination of the Na D 
abs orption (see Fig.pjb)) , very common in superwind galax- 
ies (|Heckman et al.ll2000l ). 



4.2 Fourier Analysis of the Residual Velocity 
Field 



We would like to go one step further in the analysis of the 
residual velocity field of NGC 4321. One of the most striking 
features of the residual velocity m ap is its great degree of 
symmetry. iKnapen et al. I (Il995al ) already noticed that the 
K-band image of the nuclear region of NGC 4321 presented 
a high degree of symmetry. They performed a decomposition 
of the symmetric and antisymmetric components of this im- 
age, finding that an impressive 95% of the flux is emitted 
in the symmetric component. If, as expected, the K-band 
image is a good tracer of the mass density and hence of 
the potential, it is reasonable also to expect a high degree 
of symmetry in the kinematical signature of such a poten- 
tial. We have therefore performed the same kind of decom- 
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position in our residual velocity field. The result of such a 
procedure is shown in Fig. [8] (a) , where the residual veloc- 
ity map is shown in colours with regions of symmetry and 
antisymmetry dominance plotted on white and black con- 
tours respectively. The three distinctive regions mentioned 
in Section l4~Tl now stand out more clearly. The central region 
is the only one that shows a relevant symmetric dominance 
(with blue-shifted velocities) while the bar and arm influence 
zones are regions of antisymmetric dominance. The rest of 
the map does not show any definite symmetry. 

It is not surprising that most of the residuals be- 
long to the antisymmetric component. As iKnapen et al.l 
( 1995b) showed, there is a clear m = 2 c omponen t in th e 
potential generated by the n uclear bar. ICanzia 3 (|l993h . 
and later ISchoenmakers et al.l (| 1997t ) in an analytical way, 
showed that a perturbation of harmonic number, m, in the 
potential shall generate m — 1 and m + 1 signatures in the 
velocity field (dominating inside and outside corotation re- 
spectively). 

The antisymmetric component is dominated by an har- 
monic m — 1 mode, while the symmetric component is con- 
tributed by a combination of even modes (mostly m — 0). 
To show this in a clearer way, we have also performed a 
Fourier analysis of the residual velocity field. The residual 
velocity field is decomposed into its Fourier components in 
rings of different radii. Each ring is therefore characterized 
by the amplitude F m (R) and phase <f)m(R) of each Fourier 
component m. The results are shown for the first four modes 
in Figure [8] (b). The amplitude of modes m > 4 is negligible 
at all radii. The three regions of definite symmetry coincide 
with the regions of high residuals already mentioned: m — 
dominance for the nuclear region, and m = 1 for the bar 
(with a very constant phase angle, and spiral arms (with a 
phase angle varying with radius) . The m = 3 mode is second 
in importance in the bar and spiral arms regions. 

In our case, if the bar m — 2 perturbation to the po- 
tential is responsible for the deviations from circu l ar rot a- 
tion, we expect, according to ISchoenmakers et alj (Il997l ) a 
strong m — 1 signature (with a secondary m — 3) in the 
residual velocity field (as we are inside the corotation radius 
and Qb < fi). This is exactly what is observed in our data 
for NGC 4321. Although a spiral structure is also a m=2 
distortion and would therefore produce m-l=l and m+l=3 
signatures in the velocity field, we strongly favour the bar 
scenario. The constancy of the phase angle of the m-l=l 
in the innermost region (see Figure |9j clearly indicates a 
perturbation with a constant PA, as one would expect from 
a bar. Moreover, the NIR photometry of NGC 4321 clearly 
shows the nuclear bar, but no trace of spiral arms in the 
mass distribution. 

The Fourier analysis therefore shows clearly that the ve- 
locity residuals are not local perturbations but global modes 
(due to the high degree of symmetry) generated by the 
m — 2 potential perturbation of the nuclear bar. The spiral 
arms are also generated by the bar potential. The generation 
of the spiral arms is most probably caused by shocks origi- 
nated by gas orbit cr owding at a resonance (see for example 
ISakamoto et alJll99S( ). 



100- A "a, 



Figure 9. Phase angle for m=l mode in the Fourier anal- 
ysis of the gas residual velocity field. The bar region 
(between dashed vertical lines) shows a constant phase 
and the spiral arms region (between dot-dashed vertical 
lines) shows a phase varying with radius. 



4.3 Stellar Rotation Curve. Application of the 
Tremaine— Weinberg Method 

We calculate the stellar rotation curve from the stellar ve- 
locity map (see Fig. [4] (a)) using the same method employed 
for the ionized gas. We fix the position of the centre, incli- 
nation, P.A. and systemic velocity to the values used for the 
ionized gas, and leave the circular velocity as the only free 
parameter in the fit. Fig. [4] (a) may give the wrong impres- 
sion of an offset between the kinematical and photometrical 
centers. This is due to an combination of facts, the most 
important being that the systemic velocity plotted in Fig. [4] 

(a) is the one calculated for the gas (which is used in the 
rotation model). The systemic velocity computed from the 
stellar velocity field would give a slightly lower value, re- 
sulting in a iso-velocity contour closer to the photometrical 
center. We have chosen to keep the determination of the 
systemic velocity from the ionized gas (with a much higher 
S/N) despite this spurious effect. 

The resulting stellar rotation curve is shown in Figure [4] 

(b) . It rises more slowly in the innermost region but shows 
a very similar behaviour to the gaseous rotation curve. The 
values of the stellar rotational velocity are always smaller 
than the gaseous ones. Qualitatively, this is to be expected 
owing to the stars' asymmetric drift. 

The pattern speed of the bar in NGC 4321 and the exis- 
tence of two nes ted bars or a single bar have been a much de- 
bated issue (see lHernandez et al.ll2005l . for a review). In any 
case, what seems to be clear is that there are several relevant 
pattern spe eds in NGC 4321 . Mos t observational determi- 
nations (see iHernandez et al.l [20051 . and references therein) 
seem to agree on a value of Q p « 20-35 km s _1 kpc -1 for 
the outer spiral arms. The discrepancies arise in the in- 
ner region and over whether the large scale bar and the 
nuclear bar are the same or different structures. Compar- 
ison of numerical simulations with ob servations have left 
two surviving alter n atives . On one side IWada et all {l998) 
and IKnapen et ail (|200Ch favour a single bar scenario in 
which the pattern speed of the bar is f2(, ~ 70 km s _1 kpc -1 . 
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Figure 10. Results of the TW method. This plot shows the re- 
sults for the unweighted data (empty diamonds) with the best fit 
(dashed line)). The symmetrized version is shown in solid trian- 
gles with the best fit as a solid line (see text for details). Both fits 
(for the unweighted and symmetrized data) provide very similar 
results with in a value of the slope of Q(,=160±70 km s — 1 kpc -1 . 




Figure ll.<y>vs>X> plot on top of an image of the 
NaD line strength. 



If the large scale bar has to rotate at such a high speed, 
this scenario seems to contradict the above-mentioned fact 
that most observational studies determine a much lower pat- 
tern speed for the outer regions (where both the large scale 
bar and spiral arms coexist). The a lternative model pro- 
posed bv lGarcla— Burillo et alJ (|l998l ) consists of two nested 
bars with different (although closely linked) pattern speeds. 
Their best solution provides fast and slow pattern speeds of 
f2/=160 km s -1 kpc -1 and fi 3 =23 km s -1 kpc -1 respec- 
tively. This alternative has the virtue of being able to drive 
gas inwards thanks to the fast mode (as in the single bar 
model) in the inner region while keeping the observed slow 
mode at large radii. 

The best observational determination of the pattern 
speeds of structures in galactic discs is pr o vided by the 
method proposed by iTremaine fc Weinberg! (| 1984h (here- 
after TW). This method does not rely on any particular 



theory or numerical simulations. It relies exclusively on the 
continuity equation. In a cartesian coordinated system (x, y) 
where the x- and y-axes are aligned with the galaxy's major 
and minor axis respectively, they showed that if one can find 
observationally a tracer of the mass density and the velocity 
of a galactic component that satisfies the continuity equa- 
tion, then the underlying pattern speed can be calculated 
from this observational information as: 



1 



I Zo Z(x)(v los (x) - v sys )dx i (y los } 



sini J™ Y,(x)[x - x c )dx 



(1) 



where is the observed intensity (assumed to be propor- 
tional to the surface mass density), vi os (x) is the line of sight 
velocity, v sys is the systemic velocity and x c is the x coordi- 
nate of the centre. The integrals are performed along slits or 
slices parallel to the major axis direction. TW also pointed 
out the convenience of using an odd weighting function in 
y by using, for example , slits offset by the same am ount on 
both sides of the centre. iMerrifield &: Kuiikenl l|l995t) refined 
the method by noting that a plot of (Vj os ) vs (X) provides 
a combined measurement of the pattern speed using several 
slits, which is insensitive to errors in the centre position and 
systemic velocity. 

As commented above, the TW method has been applied 
to sev eral gaseous compo nents, such as Ha ([H ernande z et al.l 
l2005h or molecular gas (|Rand fc Wallinl 12004 ). However, it 
is arguable whether these components satisfy the continu- 
ity equation. Stars, particularly old stars, on the contrary, 
are assumed to survive several cycles around the centre and 
thus, are expected to satisfy the continuity equation. It is, 
therefore, interesting to use our stellar velocity field and the 
strength of a stellar absorption line (as a tracer of the mass 
surface density) to apply the TW method to the inner region 
of NGC 4321. The integrals in Eq.[T]are performed along 21 
slits parallel to the major axis and centred on the galaxy cen- 
tre. The slits are offset by 1". The resulting (Vi 03 ) vs (X) plot 
is shown in Figure 1101 The plot shows also a symmetrized 
version by using an odd S(y — y c ) —S(y + y c ) weighting func- 
tion in the integrals. The < Y > vs < X > plot shown in 
Figure ITTI removes the possible degeneracies b etween differ- 
ent pattern speeds (see iHernandez et al.ll2005l ). 

Both fits (for the unweighted and symmetrized data) 
provide very similar results with a slope value of Qt = 
160T70 km s -1 kpc -1 . The quoted error contains the con- 
tribution from uncertainties in the P.A. As a cross-check 
w e have also used th e K s infrared image of NGC 4321 
of iKnapen et all (|2003l ) as the mass tracer for the old stars. 
The results in this case provide a somewhat higher value for 
Qb although compatible with the previous one within the 
large error bars. 

This result is very different from the one associated 
to the single bar scenario (fij, ~ 70 km s _1 kpc -1 ). How- 
ever, in spite of the numerical coincidence the agreement 
wi th the fast mode of t he ne sted bar scenario proposed 
bv lGarcfa— Burillo et al.l |l998), is not conclusive. 

A linear trend of the (Vj os ) vs (X) relationship like the 
one observed in Figure [10] is expected for the bar rotation 
pattern but it is not exclusive of it. In fact, the stellar ro- 
tation curve of the inner region of NGC 4321 is approxi- 
mately described by a linear relationship of about 100±20 
km s -1 kpc -1 between 3 and 12 arcsecs. Thus it would be 
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not surprising if the convolution (see Eq.[l]) of the corre- 
sponding solid-body like rotating velocity field with a bar of 
finite width will produce a (Vi 03 ) vs (X) linear trend with 
Q ro t ~ 100 km s _1 kpc -1 . A direct test repeating the com- 
plete procedure to apply the T-W method but now to a syn- 
thetic velocity field generated with the stellar rotation curve 
of Figure [4] yields a value Q ro t = 130 ± 30 km s _1 kpc -1 . 
Thus, Q ro t and Qb = 160±70 km s _1 kpc -1 are compatible 
within errors. In fact, owing to the relative closeness be- 
tween these two magnitudes (if the two nested bars scenario 
is correct) a robust determination of Qb would need a strong 
diminution of the errors. 

This new determination of the pattern speed has obvi- 
ously an impact on the location of the important resonances. 
A higher pattern speed moves most resonances to smaller 
radii. We have used the linear epyciclic approximation to 
make an estimation of the resonance location. There are two 
Inner Lindblad Resonances located at radii of approximately 
6 and 12 arcsec respectively. The Ultra Harmonic Resonance 
(Q — ft/4 = Qb) is located at a radius of about 14 arcsec. 
The star forming ring seen in Ha is located right inside the 
Outer Inner Lindbald Resonance and Ultra Harmonic Res- 
onance. Whether this star forming ring is related to some 
of these resonances is indeed a possibility, but we cannot be 
conclusive in this respect. In any case, these values must be 
considered very cautiously, as they depend critically on the 
shape of the rotation curve in its innermost region, and on 
the validity of the linear epyciclic approximation, which is 
probably a poor one in these regions. 

A few words of caution about the determination of the 
bar pattern speed must be said at this point. The original 
TW method eliminates some undesired integrals of spatial 
derivatives by using the fact that the integrals extend to in- 
finity, where the density goes to zero. We obviously cannot 
fulfil this requirement, as most observations having a lim- 
ited field of view can not. It is usually assumed that if the 
intensity is low enough at the borders of the field of view, 
the method is still usable. As we have seen before, there 
seems to be evidence that in NGC 4321 there are several 
speed patterns that are relevant at different radii. If this is 
the case, integrating along a slit that runs over the two pat- 
tern speed regions may contaminate the results of the TW 
method (for the inner, faster mode), even if we select the ap- 
propriate region along the y-axis. In fact, the TW method 
can be gener alized for a continuou sly varying pattern speed 
with radius l|Merrifield et alj|2006h although its application 
requires the solution of a Volterra integral equation, which 
is much more complex and needs very high quality data over 
the whole region of interest. It would then be necessary to 
have high S/N data of the type presented here covering the 
whole disc of NGC 4321 to perform such an analysis and to 
solve this issue definitively. 



5 SUMMARY 

We have performed a Fourier analysis of the residual velocity 
field of the ionized gas and determined the stellar bar pattern 
speed to study the departures from regular rotation and how 
they relate to the gravitational potential. The main results 
of this work can be summarized in the following points: 

(i) The ionized gas kinematics presents, on top of the 



galactic rotation, three clearly distinctive regions: a) A small 
nuclear region of blue-shifted velocities, b) A bar-dominated 
region where the residuals originate from the gas orbits 
aligned with the bar major axis, c) A region dominated by 
the streaming motion across two symmetric spiral arms. 

(ii) The analysis of the gas kinematical symmetries shows 
these three regions in a very clear way as regions of dom- 
inance of symmetry or antisymmetry, or, alternatively, of 
m = and m — 1 (secondary m = 3) Fourier modes. This 
is exactly what theory predicts if they are signatures of an 
m — 2 perturbation of the potential. Thus, the outstanding 
morphological and kinematic perturbations of the circumnu- 
clear region of NGC 4321 present the required (high) degree 
of symmetry to be explained by the m = 2 potential of a 
bar lying in this region. 

(iii) We have detected this (otherwise hidden) bar at op- 
tical wavelengths by means of the 2D distribution of the 
equivalent width of the stellar absorption lines correspond- 
ing to the old stellar population. This optical bar matches 
well the inner bar found in the near-infrared by Knapen et 
al. (2003). 

(iv) We apply the Tremaine-Weinberg method to the 
stellar data to estimate the pattern speed of the inner bar. 
The obtained value of Qb = 160 ± 70 km s _1 kpc -1 is very 
different from the value proposed in the simple bar scenario. 
However, the uncertainties in the p attern speed determi- 
nation prevent the confirmation of iGarcfa— Burillo et al.l 
l|l99St ) results based in the existence of two nested bars. 

(v) The nuclear blue-shifted region is probably the signa- 
ture of a nuclear outflow. This is supported by the blueward 
asymmetry of the emission line profiles and by the presence 
of a blue-shifted interstellar gas component in the Na D ab- 
sorption. 

(vi) Finally, while star formation seems to take place in 
kinematically quiet regions, the dust location coincides with 
high velocity gradients. 
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